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Abstract — In  an  ad  hoc  cognitive  radio  network,  secondary 
users  access  channels  temporarily  unused  by  primary  users  and 
the  existence  of  a  communication  link  between  two  secondary 
users  depends  on  the  transmitting  and  receiving  activities  of 
nearby  primary  users.  Using  theories  and  techniques  from  con¬ 
tinuum  percolation  and  ergodicity,  we  analytically  characterize 
the  connectivity  of  the  secondary  network  defined  in  terms  of 
the  almost  sure  finiteness  of  the  multihop  delay,  and  show  the 
occurrence  of  a  phase  transition  phenomenon  while  studying 
the  impact  of  the  temporal  dynamics  of  the  primary  traffic 
on  the  connectivity  of  the  secondary  network.  Specifically,  as 
long  as  the  primary  traffic  has  some  temporal  dynamics  caused 
by  either  mobility  and/or  changes  in  traffic  load  and  pattern, 
the  connectivity  of  the  secondary  network  depends  solely  on  its 
own  density  and  is  independent  of  the  primary  traffic;  otherwise 
the  connectivity  of  the  secondary  network  requires  putting  a 
density-dependent  cap  on  the  primary  traffic  load.  We  show  that 
the  scaling  behavior  of  the  multihop  delay  depends  critically 
on  whether  or  not  the  secondary  network  is  instantaneously 
connected.  In  particular,  we  establish  the  scaling  law  of  the 
minimum  multihop  delay  with  respect  to  the  source-destination 
distance  when  the  propagation  delay  is  negligible. 

Index  Terms — Ad  hoc  cognitive  radio  network,  connectivity, 
traffic  dynamics,  multihop  delay,  continuum  percolation,  ergod¬ 
icity. 

I.  Introduction 

In  spectrum  overlay  networks,  primary  and  secondary  users 
share  a  common  spectrum  in  a  hierarchical  manner  to  achieve 
spectrum  efficiency  and  interoperability  [1],  By  sensing  and 
learning  the  communication  environment  via  their  cognitive 
radios  [2],  secondary  users  identify  and  exploit  instantaneous 
and  local  spectrum  opportunities  while  avoiding  unacceptable 
interference  to  primary  users  [1], 

We  analytically  characterize  the  connectivity  and  multihop 
delay  of  the  secondary  network.  The  existence  of  a  commu¬ 
nication  link  between  two  secondary  users  depends  on  not 
only  their  separation  but  also  the  occurrence  of  the  spectrum 
opportunity  determined  by  the  transmitting  and  receiving 
activities  of  nearby  primary  users.  It  is  this  interaction  with  the 
primary  network  that  makes  the  problem  fundamentally  dif¬ 
ferent  from,  and  the  analysis  considerably  more  complex  than 
their  counterparts  in  homogeneous  networks.  A  qualitative  and 
quantitative  characterization  of  the  impact  of  primary  traffic  on 
the  secondary  network  is  thus  critical  for  understanding  the 
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performance  limit  of  ad  hoc  cognitive  radio  (CR)  networks, 
and  is  the  main  topic  of  this  paper. 

A.  Main  Results 

We  consider  a  Poisson  distributed  secondary  network  over¬ 
laid  with  a  Poisson  distributed  primary  network  in  an  infinite 
two-dimensional  Euclidean  space1.  We  define  connectivity  via 
the  finiteness  of  the  minimum  multihop  delay  (MMD)  between 
two  randomly  chosen  secondary  users,  referred  to  as  finite- 
delay  connectivity  (fd-connectivity).  Specifically,  the  network 
is  fd-disconnected  if  the  MMD  between  two  randomly  chosen 
secondary  users  is  infinite  almost  surely  (a.s.),  and  is  fd- 
connected  if  the  MMD  is  finite  with  a  positive  probability 
(wpp.).  Notice  that  the  MMD  considered  here  is  not  the 
multihop  delay  for  a  specific  routing  protocol.  Instead,  it  is  the 
minimum  multihop  delay  that  can  be  achieved  by  any  routing 
protocol.  The  MMD  thus  specifies  a  fundamental  performance 
limit  and  provides  a  benchmark  for  comparison. 

We  consider  temporal  dynamics  in  the  primary  traffic  which 
could  be  caused  by  mobility  and/or  changes  in  the  traffic 
load  and  pattern.  We  assume  that  the  secondary  network  is 
static.  Under  the  Poisson  model,  the  two  key  parameters  that 
characterize  the  topological  structure  of  the  secondary  network 
and  the  primary  traffic  load  are  the  density  As  of  the  secondary 
users  and  the  sequence  {A pr(t)  :  t  >  0}  of  the  densities  of 
the  primary  transmitters.  The  fd-connectivity  of  the  secondary 
network  can  thus  be  characterized  by  a  partition  of  the  infinite- 
dimensional  space  (As,{Apt(/)  :  t  >  0}). 

Although  the  above  partition  appears  to  be  intractable,  we 
show  that  as  long  as  the  primary  traffic  has  some  temporal 
dynamics  (no  matter  how  small  the  range  of  the  dynamics 
is),  the  fd-connectivity  of  the  secondary  network  depends 
solely  on  its  own  density  A s  and  is  independent  of  the 
densities  {App(i)}  of  the  primary  transmitters,  as  illustrated 
in  Fig.  1(a).  In  other  words,  no  matter  how  heavy  the  primary 
traffic  is,  the  secondary  network  is  fd-connected,  as  long  as  its 
density  A s  exceeds  the  critical  density  Ac  of  a  homogeneous 
network  (i.e.,  in  the  absence  of  the  primary  network).  Note 
that  when  A s  >  A c,  there  is  a.s.  a  unique  infinite  connected 
component  (ICC)  [3,  Chapter  3]  in  the  secondary  network 
formed  by  topological  links  (a  topological  link  exists  between 
two  users  that  are  within  communication  range).  We  show  that 
for  any  two  secondary  users  in  this  ICC,  the  MMD  is  finite 
a.s.  The  intuition  is  that  messages  can  traverse  a  topological 
path  connecting  the  two  secondary  users  by  making  stops 

'This  infinite  network  model  is  equivalent  in  distribution  to  the  limit  of  a 
sequence  of  finite  networks  with  a  fixed  density  as  the  area  of  the  network 
increases  to  infinity,  i.e.,  the  so-called  extended  network. 
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Fig.  1.  FD-Connectivity  of  ad  hoc  CR  networks  when  the  primary  traffic  has  temporal  dynamics  (a)  and  no  temporal  dynamics  (b).  The  critical  density 
AJ  of  the  secondary  users  is  defined  as  the  infimum  density  of  the  secondary  users  that  ensures  instantaneous  connectivity  under  a  positive  density  of  the 
primary  transmitters,  and  is  equal  to  the  critical  density  Ac  of  a  homogeneous  network;  the  upper  boundary  ApT(A s)  is  defined  as  the  supremum  density  of 
the  primary  transmitters  that  ensures  instantaneous  connectivity  with  a  fixed  density  As  of  the  secondary  users. 


in  between  to  wait  for  spectrum  opportunities,  and  more 
importantly,  the  waiting  time  is  finite  a.s.  due  to  the  temporal 
dynamics  of  the  primary  traffic.  Since  the  percentage  of  the 
secondary  users  in  this  ICC  is  strictly  positive,  it  follows  that 
the  MMD  between  two  randomly  chosen  secondary  users  is 
finite  wpp.,  i.e.,  the  secondary  network  is  fd-connected. 

On  the  other  hand,  when  the  primary  network  is  static,  we 
show  that  the  secondary  network  is  fd-connected  if  and  only 
if  (iff.)  it  is  instantaneously  connected,  as  shown  in  Fig.  1(b). 
The  secondary  network  is  instantaneously  connected  if  it  has  a 
unique  ICC  formed  by  communication  links  a.s.  The  existence 
of  a  communication  link  requires  the  existence  of  a  topological 
link  and  the  presence  of  a  spectrum  opportunity  determined 
by  the  transmitting  and  receiving  activities  of  nearby  primary 
users.  Due  to  this  requirement,  the  instantaneous  connectivity 
puts  a  cap  on  the  tolerable  primary  traffic  which  is  an 
increasing  function  of  the  density  A s  of  the  secondary  users 
(see  Fig.  1(b)).  Moreover,  given  a  static  primary  network, 
the  set  of  communication  links  in  the  secondary  network  is 
fixed  over  time.  It  implies  that  if  a  topological  link  does  not 
see  an  opportunity  at  the  beginning,  then  it  will  never  see 
it.  Thus,  messages  from  one  secondary  user  can  only  reach 
another  secondary  user  within  the  same  connected  component 
formed  by  communication  links.  If  the  secondary  network  is 
instantaneously  connected,  then  wpp.  two  randomly  chosen 
secondary  users  belong  to  this  ICC  formed  by  communication 
links2,  and  the  MMD  between  them  is  finite;  otherwise  they 
belong  to  two  different  finite  connected  components  a.s.,  and 
they  are  inaccessible  from  each  other,  i.e.,  the  MMD  is  infinite. 

Although  the  primary  traffic  does  not  affect  the  fd- 
connectivity  of  the  secondary  network  when  it  has  temporal 
dynamics,  it  does  affect  the  behavior  of  the  MMD.  Indeed,  we 
show  that  the  scaling  behavior  of  the  MDD  with  respect  to  the 
source-destination  distance  is  starkly  different  depending  on 
whether  the  secondary  network  is  instantaneously  connected 
wpp.  or  not.  Notice  that  the  multihop  delay  in  the  secondary 
network  consists  of  two  components:  the  propagation  delay 
and  the  waiting  time  at  each  hop  for  the  occurrence  of  a 
spectrum  opportunity.  When  the  propagation  delay  is  negligi- 

-  [t  is  shown  in  [4]  that  there  exists  either  zero  or  one  ICC  formed  by 
communication  links  in  the  secondary  network  a.s. 


ble,  we  show  that  if  the  secondary  network  is  instantaneously 
connected  wpp.,  the  MMD  is  asymptotically  independent  of 
the  source-destination  distance;  otherwise  the  MMD  scales  at 
least  linearly  with  the  source-destination  distance.  We  also 
study  the  case  of  nonnegligible  propagation  delay.  Simulations 
show  that  the  MMD-to-distance  ratio  for  a  secondary  network 
that  is  instantaneously  connected  wpp.  can  be  orders  of 
magnitude  smaller  than  that  for  a  secondary  network  that  is 
not  instantaneously  connected  a.s. 

These  analytical  results  also  provide  important  insights 
and  design  guidelines  for  practical  systems.  Since  almost  all 
primary  networks  have  temporally  dynamic  traffic,  it  follows 
from  the  result  on  fd-connectivity  that  the  accessibility  be¬ 
tween  two  secondary  users  is  independent  of  the  presence  of 
the  primary  network,  although  it  may  incur  a  larger  multihop 
delay.  From  the  result  on  multihop  delay,  we  can  see  that 
if  the  primary  network  has  heavy  traffic,  then  the  secondary 
network  can  only  be  used  for  delay-tolerant  applications; 
conversely,  if  a  secondary  network  is  deployed  for  delay- 
sensitive  applications,  then  it  should  be  operated  within  the 
instantaneous  connectivity  region  Cins  for  a  positive  portion 
of  time3,  which  imposes  restrictions  on  the  traffic  load  of  the 
primary  network  or  on  the  density  of  the  secondary  network. 

B.  Related  Work 

There  have  been  only  a  few  results  on  the  connectivity 
of  ad  hoc  CR  networks.  The  Laplacian  matrix  is  used  to 
approximately  characterize  the  graph  connectivity  in  [5];  but 
this  does  not  characterize  the  multihop  delay,  and  it  does  not 
take  into  account  the  impact  of  the  receiving  activities  of  the 
primary  network  on  the  secondary  network. 

Different  types  of  connectivity  of  homogeneous  networks 
{i.e.,  secondary  network  only)  have  been  well  studied  in  [6-14] 
and  references  therein.  The  theory  of  continuum  percolation 
has  been  used  by  Dousse  et  al.  in  analyzing  the  connectivity 
under  the  worst  case  mutual  interference  [10,  11].  In  [12,  13], 
the  connectivity  and  the  multihop  delay  in  a  homogeneous 

3 Since  {A pt{1)}  is  ergodic,  it  follows  from  Fact  Al  (in  Appendix  A)  that 
the  instantaneous  connectivity  of  the  secondary  network  wpp.  is  equivalent  to 
that  of  the  secondary  network  for  a  positive  portion  of  time. 
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network  with  static  or  dynamic  on-off  links  are  investigated 
from  a  percolation-based  perspective.  A  dynamic  connectivity 
graph  for  ALOHA  networks  is  introduced  in  [14]  to  establish 
the  scaling  law  of  the  delay  with  respect  to  the  source- 
destination  distance. 

In  [15],  we  have  studied  the  connectivity  and  multihop 
delay  of  ad  hoc  CR  networks  under  the  assumption  that  the 
realizations  of  the  primary  network  are  i.i.d.  across  slots,  but 
we  had  not  obtained  the  necessary  and  sufficient  condition 
for  the  independence  of  the  fd-connectivity  of  the  secondary 
network  from  the  primary  traffic.  In  this  paper,  not  only  is 
this  i.i.d.  assumption  replaced  by  a  more  realistic  assumption 
under  which  the  realizations  of  the  primary  network  can  be 
temporally  correlated,  but  also  the  necessary  and  sufficient 
condition  is  provided.  As  detailed  in  Sec.  III-D  and  Sec.  IV, 
the  relaxation  of  this  i.i.d.  assumption  significantly  complicates 
the  analysis,  especially  the  one  for  fd-connectivity. 

II.  Network  Model 

We  consider  a  Poisson  distributed  secondary  network  over¬ 
laid  with  a  Poisson  distributed  primary  network  in  an  infinite 
two  dimensional  Euclidean  space.  The  primary  network  adopts 
a  synchronized  slotted  structure  with  slot  length  Tg.  Thus  Tg 
can  be  considered  to  be  the  time  constant  of  the  spectrum 
opportunities  which  are  determined  by  the  transmitting  and 
receiving  activities  of  the  primary  users.  Without  loss  of 
generality,  we  set  Tg  =  1. 

At  t  =  0,  primary  transmitters  are  distributed  according 
to  a  two-dimensional  Poisson  point  process  npr(0)  with 
density  Apt(0).  At  each  t  >  0,  each  primary  transmitter  has 
a  probability  q(t)  of  not  transmitting  at  t  +  1,  due  to  lack  of 
packets  or  perceived  channel  conditions.  The  primary  transmit¬ 
ters  which  continue  transmitting  will  move  to  a  new  position 
at  t  +  1  according  to  a  random  displacement  vector  m(t ) 
with  a  finite  variance  in  each  direction.  The  movements  are 
i.i.d.  for  different  primary  transmitters,  and  for  each  primary 
transmitter,  m(t)  is  i.i.d.  across  slots4.  Based  on  the  thinning 
theorem  and  [17] [Proposition  1.3],  the  primary  transmitters 
which  continue  transmitting  will  form  a  two-dimensional 
Poisson  point  process  with  density  App(f)[l  —  g(f)]  at  t  +  1. 

Let  nPT(f)  (t  >  0)  denote  the  point  process  of  the 
primary  transmitters  at  t.  At  t  +  1,  some  primary  nodes, 
that  were  silent  in  slot  t,  may  start  transmitting.  They  are 
distributed,  independent  of  n pp(f),  according  to  a  two- 
dimensional  Poisson  point  process  with  density  p(f)App(f), 
where  the  multiplicative  factor  p(t)  >  0  may  exceed  unity5. 
It  follows  by  induction  that  the  point  process  of  transmitters 
IIp7’(f  +  1)  containing  the  old  primary  transmitters  at  t  and 
the  new  primary  transmitters  at  t  +  1  is  Poisson  with  density 
Xpr(t  +  1)  =  App(f)[l  —  q(t)  +  p(t)].  The  random  process 


{App(f)}  is  assumed  to  be  stationary  and  ergodic.  The  two 
related  random  processes  {q(t)  £  [0, 1]}  and  {p(t)  >  0}  are 
assumed  to  be  stationary  and  ergodic;  they  may  be  correlated 
with  {A P7’(f)}. 

The  primary  receivers  are  randomly  (may  not  be  uniformly) 
located  within  the  transmission  range6  Rp  of  their  corre¬ 
sponding  transmitters  at  each  t,  and  their  relative  positions 
with  respect  to  their  corresponding  transmitters  can  be  either 
fixed  or  a  stationary  and  ergodic  random  process  over  time. 
Based  on  the  displacement  theorem  [18,  Chapter  5],  it  can 
be  shown  that  at  each  t,  the  primary  receivers  form  another 
two-dimensional  Poisson  point  process  n pp(t)  with  density 
A pp(f),  which  is  correlated  with  n prit). 

Secondary  users  are  distributed  according  to  a  two- 
dimensional  Poisson  point  process  ng  with  density  A g,  which 
is  independent  of  {npp(f)}  and  {npp(t)}.  The  locations 
of  the  secondary  users  are  static  over  time,  and  they  have 
a  uniform  transmission  range  rp. 

III.  Connectivity 

In  this  section,  we  analytically  characterize  the  connectivity 
of  the  secondary  network.  In  particular,  we  show  the  occur¬ 
rence  of  a  phase  transition  phenomenon  in  terms  of  the  impact 
of  the  temporal  dynamics  of  the  primary  traffic  on  the  fd- 
connectivity  of  the  secondary  network. 

A.  Topological  Link  vs.  Communication  Link 

A  topological  link  exists  between  any  two  secondary  users 
that  are  within  each  other’s  transmission  range.  Thus,  topolog¬ 
ical  links  in  the  secondary  network  are  independent  of  the  pri¬ 
mary  network.  In  contrast,  as  discussed  in  the  next  paragraph, 
the  existence  of  a  communication  link  between  two  secondary 
users  depends  not  only  on  the  distance  between  them  but 
also  on  the  availability  of  the  communication  channel,  i.e.,  the 
presence  of  a  spectrum  opportunity.  As  a  result,  even  in  a  static 
secondary  network,  communication  links  are  time-varying  due 
to  the  temporal  dynamics  of  spectrum  opportunities. 


Fig.  2. 
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4This  assumption  of  “i.i.d.  across  slots”  is  made  so  that  we  can  use  the 
classical  central  limit  theorem  (see  Appendix  B).  Since  the  central  limit 
theorem  can  be  extended  to  the  two  cases  of  identical  but  weakly  dependent 
distributions  and  martingales  [16,  Sec.  7.7],  this  assumption  can  be  relaxed. 

5  We  introduce  p(t)  for  convenience.  One  can  consider  this  as  a  birth-death 

process;  nodes  die  when  they  have  no  more  packets  to  send,  and  are  (re-)bom 
when  they  do.  We  could  also  consider  this  from  a  duty-cy cling  perspective: 
nodes  sleep  and  wake  up. 


We  consider  the  disk  signal  propagation  and  interference 
model  as  illustrated  in  Fig.  2.  There  exists  an  opportunity  from 
/i,  the  secondary  transmitter,  to  is,  the  secondary  receiver,  if  the 
transmission  from  p  does  not  interfere  with  primary  receivers 

6Here  we  assume  that  all  the  primary  transmitters  use  the  same  power  and 
the  transmitted  signals  undergo  isotropic  path  loss. 
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in  the  solid  circle,  and  the  reception  at  v  is  not  affected  by 
primary  transmitters  in  the  dashed  circle  [19].  Referred  to  as 
the  interference  range  of  secondary  users,  the  radius  77  of  the 
solid  circle  centered  at  p  depends  on  the  transmission  power  of 
p  and  the  interference  tolerance  of  primary  receivers,  whereas 
the  radius  Rj  of  the  dashed  circle  (the  interference  range  of 
primary  users)  depends  on  the  transmission  power  of  primary 
users  and  the  interference  tolerance  of  the  secondary  user  v. 

It  follows  from  the  above  discussion  that  spectrum  op¬ 
portunities  are  asymmetric.  Specifically,  a  channel  that  is  an 
opportunity  when  p  is  the  transmitter  and  v  the  receiver  may 
not  be  an  opportunity  when  v  is  the  transmitter  and  p  the 
receiver.  Since  unidirectional  links  are  difficult  to  utilize,  es¬ 
pecially  for  applications  with  guaranteed  delivery  that  require 
acknowledgements,  we  only  consider  bidirectional  links  in  the 
secondary  network  when  we  define  connectivity. 


B.  Instantaneous  Connectivity  vs.  Topological  Connectivity 

In  each  slot  t,  we  can  obtain  an  undirected  random  graph 
Qh (As,  Xprit),  t)  consisting  of  all  the  secondary  users  and 
their  communication  links,  which  represents  the  instantaneous 
connectivity  of  the  secondary  network  in  this  slot.  As  illus¬ 
trated  in  Fig.  3,  this  graph  Gh{Xs,  Xpr{t),t)  is  determined 
by  the  three  Poisson  point  processes  in  slot  t, :  II5,  Ylprit), 
and  Ilpp(f),  where  II pr{t)  and  Ilpp(f)  are  correlated. 


Fig.  3.  A  realization  of  the  random  graph  Gh{^S,  A pxif),  t)  which  consists 
of  all  the  secondary  users  and  their  communication  links  in  slot  t  (denoted  by 
solid  lines).  The  solid  circles  denote  the  interference  regions  of  the  primary 
transmitters  within  which  secondary  users  cannot  successfully  receive,  and  the 
dashed  circles  denote  the  required  protection  regions  for  the  primary  receivers 
within  which  secondary  users  should  refrain  from  transmitting. 


We  define  the  instantaneous  connectivity  of  the  secondary 
network  in  slot  t  as  the  a.s.  existence  of  a  unique  ICC 
in  Gh[Xs,  App(f),  t).  Given  the  transmission  power  and  the 
interference  tolerance  of  both  the  primary  and  the  secondary 
users  (i.e.,  Rp,  Ri,  rp,  and  77  are  fixed),  the  instantaneous 
connectivity  region  Cins  f°r  slot  t  is  defined  as 

C/jvs  =  {(As,  Xprit)).  Gh (As,  Xprit),  t)  is  connected}.  (1) 
The  upper  boundary  \*PT(\ g)  of  Cins  is  defined  as 

App(As)  =  sup{Apr(f)  :  f?(As,  Apt (t),  t)  is  connected.}.  (2) 
The  critical  density  of  the  secondary  users,  X*s  is  defined  as 

A *s  =  inf{As  :  3App(f)  >  0  s.t.  Gh{Xs,  Xprit), t)  is  connected}. 


It  is  shown  in  [4]  that  Ap  equals  the  critical  density  Ac  of  a 
homogeneous  network.  A  detailed  analytical  characterization 
of  Cins  is  given  in  [4].  Let  0(As,  Apt(/))  denote  the  prob¬ 
ability  that  an  arbitrary  secondary  user  belongs  to  the  ICC  in 
Gh{Xs,  Xprit), t),  if  one  exists,  then  we  have  that 


0(As,  App(f)) 


>  0,  if  (As,  Xprit))  e  Cins', 
=  0,  otherwise. 


(3) 


The  fd-connectivity  of  the  secondary  network  is  defined 
by  the  finiteness  of  the  MMD  between  two  randomly  chosen 
secondary  users.  To  ensure  finiteness  of  the  multihop  delay 
between  two  secondary  users,  it  is  necessary  to  have  a  path 
formed  by  topological  links  between  them,  otherwise  they  are 
not  accessible  from  each  other.  Consider  an  undirected  random 
graph  Gsi^s)  consisting  of  all  the  secondary  users  and  their 
topological  links.  Notice  that  Gsi^s)  depends  only  on  the 
Poisson  point  process  Fig  of  the  secondary  network.  Define 
the  topological  connectivity  of  the  secondary  network  as  the 
a.s.  existence  of  a  unique  ICC  in  Gsi^s)-  It  follows  that  fd- 
connectivity  implies  topological  connectivity,  i.e.,  topological 
connectivity  is  usually  weaker  than  fd-connectivity.  On  the 
other  hand,  it  is  easy  to  show  that  instantaneous  connectivity 
is  usually  stronger  than  fd-connectivity. 


C.  Connectivity  with  Static  Primary  Network 

Consider  a  static  primary  network,  i.e.,  the  sets  of  the 
primary  transmitters  and  receivers  do  not  change  over  time, 
and  their  positions  are  also  fixed.  Then,  as  shown  below, 
the  necessary  and  sufficient  condition  for  the  connectivity 
of  the  secondary  network  is  its  instantaneous  connectivity 
(see  Fig.  1(b)  for  an  illustration),  i.e.,  the  connectivity  of 

Gh(Xs,  Xpr(t),  t). 

Proposition  1:  Given  a  static  primary  network,  i.e.,  p(t)  = 
q(t)  =  0  and  m(t)  =  0  V  t,  a  necessary  and  sufficient 
condition  for  the  connectivity  of  the  secondary  network  is 
given  by  (A 5,  A_pt)  E  Cins  (defined  by  (1)),  i.e.,  the  MMD 
is  finite  wpp.  iff.  the  network  is  instantaneously  connected. 

Proof:  If  (As,  Apt)  E  Cins,  then  there  exists  a  unique 
ICC  formed  by  communication  links  a.s  in  the  secondary 
network.  It  implies  that  two  randomly  chosen  secondary  users 
belong  to  this  ICC  wpp.  Since  the  set  of  communication  links 
do  not  change,  it  follows  that  the  MMD  between  them  is  finite 
wpp.,  i.e.,  the  secondary  network  is  connected. 

If  (As,  Apt)  ^  Cins,  then  only  finite  connected  compo¬ 
nents  formed  by  communication  links  exist  a.s.  Thus,  two 
randomly  chosen  secondary  users  belong  to  two  different 
connected  components  a.s.,  and  the  MMD  between  them  is 
infinite  a.s.  ■ 


D.  Connectivity  with  Dynamic  Primary  Network 

Let  r(t)  denote  the  magnitude  of  the  displacement  vec¬ 
tor  m(f).  Consider  a  dynamic  primary  network,  where  the 
dynamics  can  be  caused  by  mobility  (E[r2(f)]  >  0)  and/or 
changes  in  traffic  load  and  pattern  (E[g(f)]  >  0  implying 
that  E[p(f)]  >  0).  As  illustrated  in  Fig.  1(a),  we  show  in  the 
following  proposition  that  a  necessary  and  sufficient  condition 
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for  the  connectivity  of  the  secondary  network  is  its  topological 
connectivity,  i.e.,  the  connectivity  of  Gs(^s)- 

Proposition  2:  Consider  a  dynamic  primary  network,  i.e., 
E[q(f)]  >  0  or  E[r2(f)]  >  0.  A  necessary  and  sufficient 
condition  for  the  connectivity  of  the  secondary  network  is 
As  >  Ac,  where  Ac  is  the  critical  density  of  homogeneous 
networks.  ■ 

Remark.  Since  the  multihop  delay  is  a  finite  sum  of  single¬ 
hop  delays,  the  a.s.  finiteness  of  the  MMD  is  implied  by  that 
of  the  single-hop  delay.  Under  the  primary  network  model 
where  its  realizations  are  correlated  across  slots,  the  single¬ 
hop  delay  is,  however,  difficult  to  analyze.  In  the  proof,  we 
use  theories  and  techniques  from  ergodic  theory  to  overcome 
this  difficulty.  Specifically,  we  establish  the  ergodicity  of  a 
measure -preserving  (m.p.)  dynamical  system  which  consists  of 
the  probability  space  associated  with  the  primary  transmitters 
and  an  m.p.  shift  transformation  in  the  time  domain.  A  brief 
introduction  to  ergodic  theory  can  be  found  in  Appendix  A. 

Proof:  If  As  <  Ac,  then  there  does  not  exist  an  infinite 
topologically  connected  component  in  Gsi^s)  a.s.  It  follows 
that  two  randomly  chosen  secondary  users  belong  to  two 
different  topologically  connected  components  a.s.,  and  the 
MMD  between  them  is  infinite  a.s. 

If  A s  >  Ac,  then  there  exists  a  unique  infinite  topologically 
connected  component  Ct  in  Gs (As)  a.s.  It  follows  that  two 
randomly  chosen  secondary  users  //  and  v  belong  to  Ct  wpp. 
In  other  words,  we  can  find  a  topological  path  L  with  finite 
hops  from  p  to  v  wpp.  Since  the  MMD  t(fi,  v)  is  bounded 
above  by  the  multihop  delay  tL{p,  v)  along  the  path  L,  it 
suffices  to  show  the  a.s.  finiteness  of  tL(p,,v),  which  is  a 
direct  consequence  of  the  following  lemma. 

Lemma  1:  Let  ts(w i,  wf)  denote  the  single-hop  delay  from 
W\  to  w 2,  where  w\  and  W2  are  connected  via  a  topological 
link.  If  E[g(f)]  >  0  or  E[r2(f)]  >  0,  then  ts{w\,W2)  <  oo  a.s. 

Proof  of  Lemma  1:  Assume  that  the  propagation  delay 
t  <  Ts  =  1  so  that  the  spectrum  opportunity  lasts  long  enough 
to  ensure  the  success  of  the  transmission.  Also  assume  that 
w\  intends  to  transmit  the  message  at  t  =  0.  Thus,  ts(wi ,  W2) 
is  the  waiting  time  tsw(wi,W2)  for  the  presence  of  the  first 
bidirectional  opportunity  plus  the  propagation  delay  r,  i.e., 

ts(wi,w2)  =  tsw  +  t  =  argmin  {IB(t)  =  1}  +  r, 
te{o,i,2,...} 

where  le(f)  is  the  indicator  of  the  event  that  a  bidirectional 
opportunity  exists  in  the  tth  slot. 

Next  we  show  the  a.s.  finiteness  of  tsw.  Let  I(u>,  d,  rx/tx) 
denote  the  event  that  there  exist  primary  receivers/transmitters 
within  distance  d  of  the  secondary  user  vu,  and  I(iu,  d,  rx/tx) 
the  complement  of  I(u;,  d,  rx/tx).  The  occurrence  of  the  bidi¬ 
rectional  opportunity  E  is  given  by 

£'  =  I(u>i,r/,rx)  n  I(tui,  Ri,  tx)  fl  I(w2,  r/,  rx)  D  I(w2,  Ri,  tx). 


Let  O  be  the  midpoint  of  the  segment  connecting  w\  and 
u>2.  Define  the  event  F  as 

F  =  I(0,  I?m,  tx), 


where  Rm  =  max  {77  +  Rp  +  ^ ,  Ri  +  ^f}.  Let  tsw  be  the 
waiting  time  for  the  first  occurrence  of  the  event  F,  i.e., 

tsw  =  argmin  (If(£)  =  1}, 
te{o,i,2,...} 

where  Ip(f)  is  the  indicator  of  the  event  F  during  the  tth  slot. 
Since  F  C  E,  we  have  tsw  <  tsw.  Thus,  we  can  show  the  a.s. 
finiteness  of  tsw  by  proving  the  a.s.  finiteness  of  tsw. 

Consider  the  stationary  random  process  {Ilpj^f)  :  f>0} 
where  Ilpp(f)  is  the  Poisson  point  process  formed  by  the 
primary  transmitters  in  slot  t.  Based  on  a  trivial  generalization 
of  the  Kolmogorov  extension  theorem,  we  can  construct  a 
double-sided  stationary  random  process  { 1 1  />p  (f)  :  f  6  Z} 
which  has  the  same  finite  dimensional  distributions  as 
{Ilpp(f)  :  t  >  0}.  Let  (flpT^PT,  Ppt)  be  the  probability 
space  of  (Ilpp(f)  :  t  £  Z}.  Let  irt  (t  £  Z)  denote  the 
realization  of  Ilpp(f).  V  cu  =  {...,  7r_i,  7T0,  7Ti,  ...}  €  GIpt, 
define  a  shift  transformation  T  as 


(Tuj)t  —  Ttt+ti  V  f  G  Z, 


(4) 


where  ( Tu>)t  denotes  the  i-th  realization  of  Tu.  Since 
{Ilpp(f)  :  f  GZ}  is  time-stationary7,  it  follows  that 
(f }pt,3pt,Ppt,T)  constitute  an  m.p.  dynamical  system. 
If  (LlpTi^pT,  Ppt,T)  is  ergodic,  which  will  be  shown  in 
Lemma  2,  it  follows  from  Fact  A1  that  a.s. 

w  n—  1  n— 1 

lim  -Vlp(fc)  =  lim  -V  Tfelp(0)  =E[Ip(0)] 

n, — >00  n  ‘  ^  n, — >00  n.  '  ^ 


k—0 


k—0 


roo 

/  exP  (— AptttRm)  dF(App)  >  0, 
Jo 


where  F(App)  is  the  CDF  of  App(0).  Thus,  tsw  <  oo  a.s.  ■ 

Now  we  only  need  to  prove  the  following  lemma  to  com¬ 
plete  the  proof  of  Lemma  1 . 

Lemma  2:  Let  (GIpt^pt,  Ppt)  be  the  probability  space 
of  {Ilpp(f)  :  t  £  Z}  and  T  the  shift  transformation  de¬ 
fined  by  (4).  If  E[g(£)]  >  0  or  E[r2(f)]  >  0,  then 

(GIpt,3pt,Ppt,T)  is  ergodic. 

Sketch  Proof  of  Lemma  2:  We  show  that 

(GlpTj^PT,  Ppt,T)  is  mixing,  which  implies  its 
ergodicity  [20,  Proposition  2.5.1].  This  is  done  by  proving 
the  asymptotic  independence  of  one  event  from  another 
transformed  event.  For  details,  please  see  Appendix  B.  ■ 
Remark.  If  E[q(t)]  =  E[p(£)]  =  0  and  E[r2(f)]  =  0, 
i.e.,  the  primary  network  is  static,  it  can  be  easily  shown 
that  the  m.p.  dynamical  system  (GIpt^pt,  Ppt,T)  is  not 
ergodic.  Consider  the  following  counterexample:  let  /  £ 
L1(GlpT13pT,  Ppt)  be  the  indicator  function  of  the  event 
that  there  does  not  exist  any  primary  transmitter  within  the 
unit  square  B\  centered  at  the  origin  at  t  =  0,  then  for  any 
u>  £  Glpp,  the  time  average  /  of  /  is  given  by 

TL  —  1  ^ 

1,  if  no  primary  tx  in  B\  at  t  =  0 
0,  otherwise. 


1  r 

f=  lim  -J2f{Tku;)  =  i 

n— >oo  ji  •4 — ' 

k= 0  k 

But  the  ensemble  average  E[/]  of  /  is  given  by 

E[/]  =  f  f{oj)dPPT  =  exp(-A  PT). 
J  17  PT 


7The  temporal  stationarity  of  {npx-(f)}  is  shown  by  [17,  Proposition  1.3]. 


6 


Based  on  Fact  Al,  we  have  that  the  m.p.  dynamical  system  in 
the  static  case  is  not  ergodic.  ■ 

Combining  Proposition  1  and  2,  we  obtain  a  necessary  and 
sufficient  condition  for  the  independence  of  the  fd-connectivity 
of  the  secondary  network  with  the  primary  network. 

Theorem  1:  Let  Ac  be  the  critical  density  of  homogeneous 
networks.  Then  A 5  >  Ac  is  a  necessary  and  sufficient  con¬ 
dition  for  the  fd-connectivity  of  the  secondary  network  iff. 
E[g(t)]  >  0  or  E[r2(f)]  >  0. 

This  theorem  implies  the  occurrence  of  a  phase  transition 
phenomenon  in  the  necessary  and  sufficient  condition  for  the 
fd-connectivity.  Specifically,  if  the  primary  network  is  static, 
the  connectivity  of  the  secondary  network  is  equivalent  to  its 
instantaneous  connectivity  which  depends  on  both  its  topology 
and  the  primary  traffic;  if  the  primary  network  is  dynamic, 
the  connectivity  of  the  secondary  network  is  equivalent  to  its 
topological  connectivity  which  depends  solely  on  its  topology 
and  is  independent  of  the  primary  traffic. 

IV.  Multihop  Delay 

In  this  section,  we  analytically  characterize  the  scaling 
behavior  of  the  MMD  with  respect  to  the  source-destination 
distance  when  the  primary  traffic  is  dynamic.  Let  C(Qs(\s)) 
be  the  ICC  in  Gs(^s)  when  A s  >  Ac  i.e.,  the  secondary 
network  is  fd-connected.  We  seek  to  establish  the  scaling  law 
of  the  MMD  between  two  arbitrary  users  in  C(Gs(^s))  with 
respect  to  the  distance  between  them.  As  shown  below,  the 
scaling  behavior  of  the  MMD  is  determined  by  whether  the 
secondary  network  is  instantaneously  connected  wpp.  or  not. 


A.  Negligible  Propagation  Delay 

When  the  propagation  delay  r  =  0,  once  a  user  has  re¬ 
ceived  the  message,  it  can  spread  the  message  instantaneously 
throughout  the  connected  component,  formed  by  communica¬ 
tion  links,  which  contains  it.  Thus,  if  the  secondary  network 
is  instantaneously  connected  during  some  time  slot,  the  source 
can  route  its  message  via  the  ICC  such  that  the  message  can 
move  a  large  number  of  hops  towards  the  destination  within 
this  slot,  leading  to  the  multi-hop  delay  being  asymptotically 
independent  of  the  source-destination  distance.  On  the  other 
hand,  if  the  secondary  network  is  always  not  instantaneously 
connected,  the  message  can  move  forward  only  a  limited 
number  of  hops  within  each  slot,  which  results  in  the  linear 
scaling  of  the  MMD.  We  state  this  formally  next. 

Theorem  2:  Assume  that  r  =  0,  and  E[g(i)]  >  0  or 
E[r2(f)]  >  0.  For  any  two  secondary  users  p,u£  C(Gs(Xs)), 
the  ICC  of  Gs (As),  let  t(p,  v)  denote  the  MMD  from  p  to  v 
and  d(p ,  v)  the  distance  between  p  and  v\  then 

T2.1  if  Pr{Apr(f)  <  ApT(Ag)}  >  0  where  \*PT(Xs)  is 
defined  in  (2), 


lim 

►  00 


g(d(p,u)) 


=  0  a.s., 


where  g(d)  is  any  monotonically  increasing  function  of 
d  with  lim  g(d)  =  00; 

d — >00 


T2.2  if  Pr{ApT(t)  <  A PT}  =  0  for  some  \'PT  >  ApT(A s), 


lim  inf 

d(n,is)—>  00 


Ht{p,v)} 

d{p,u) 


>  0. 


Pr{A PT(t)  <  ApT(As)}  >  0  implies  Pr{(As,  App(f))  € 
Ci ns}  >  0,  and  Pr{App(f)  <  XPT}  =  0  for  some  X'PT  > 
App(A s)  implies  Pr{(As,  App(f))  £  Cins}  =  0,  but  not  vice 
versa.  We  state  the  two  conditions  in  the  above  way,  because 
we  have  not  been  able  to  establish  whether  the  boundary  point 
(As,  ApT(As))  £  Cins- 

Proof  Sketch:  For  T2.1,  we  use  the  ICC  in 

Gh{Xs,  XpT(fo),to)  during  some  slot  1 0  to  construct  a  path 
from  p  to  v  such  that  the  multihop  delay  along  this  path 
is  independent  of  the  distance  d(p,  v)  (see  Fig.  4  for  an 
illustration).  Let  to  be  the  first  slot  such  that  p  belongs  to  the 
ICC  C(to)  of  Gh(Xs,  XpT(to),to),  and  wu  the  user  in  C(to) 
which  is  closest  to  v.  Since  the  propagation  delay  r  =  0,  the 
multihop  delay  from  p  to  wv  is  zero.  It  follows  that 


t{p,v)  =  to  +t{wv,v). 

Then  it  suffices  to  show  that  to  and  t(wv,  v)  are  independent 
of  d(p,  v),  which  we  prove  by  using  continuum  percolation 
theory  and  ergodic  theory.  For  details,  please  see  Appendix  C. 


Fig.  4.  An  illustration  of  the  constructed  path  Lp  from  ji  to  u  when 
Pr{(As,  Apr(t))  £  CINs}  >  0.  C(to)  is  the  ICC  of  6(As,  App(fo)) 
which  first  contains  //,  and  wv  is  the  user  in  C(fo)  which  is  closest  to  v. 


T2.2  is  proven  by  using  a  coupling  argument  [3,  Chapter 
2]  and  then  deriving  a  lower  bound  on  by  considering 

the  fact  that  the  message  from  p  can  traverse  only  a  finite 
distance  towards  v  during  each  slot.  For  details,  please  see 
Appendix  D.  ■ 

From  the  above,  we  see  that  the  existence  of  the  giant 
connected  component  can  significantly  reduce  the  multihop 
delay,  especially  when  the  destination  is  far  away  from  the 
source. 


B.  Nonnegligible  Propagation  Delay 

When  the  propagation  delay  r  >  0,  it  takes  at  least  time 
r  for  the  message  to  traverse  a  distance  rp,  which  imposes 
a  lower  bound  t  j rp  on  the  ratio  of  the  MMD  to  the  source- 
destination  distance.  This  implies  that  the  MMD  scales  at  least 
linearly  with  the  source-destination  distance. 

The  positive  propagation  delay  r  also  imposes  an  upper 
bound  Ts/t  on  the  maximum  number  of  hops  that  the 
message  can  traverse  in  a  slot  Ts-  If  the  secondary  network  is 
instantaneously  connected  in  this  slot,  this  upper  bound  can  be 
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(c)  instantaneously  connected  (\pt  =  5km  2,  r  =  0.01s)  (d)  not  instantaneously  connected  (Apr  =  30km  2,  r  =  0.01s) 

Fig.  5.  MMD-to-distance  ratio  (in  logarithmic  scale)  vs.  source-destination  distance  for  random  walk  model.  Notice  that  the  MMD-to-distance  ratio  is 
obtained  in  one  Monte  Carlo  run.  The  secondary  users  are  distributed  within  a  square  [—5km,  5km]  X  [—5km,  5km]  with  density  A g  =  700km-2.  Given  the 
transmission  range  rp  =  50 m  of  the  secondary  users.  Ay  is  larger  than  the  critical  density  Ac(50)  =  576km-2.  Other  simulation  parameters  are:  rj  =  80m, 
Rp  =  50m,  Rj  =  80m,  Tg  =  Is,  p^w  =  0.05,  rm  =  5m,  tm  =  30m. 


actually  attained  in  the  ICC  consisting  of  communication  links. 
Otherwise,  this  upper  bound  may  not  be  attained  due  to  the 
limited  diameter  of  the  finite  connected  components  formed  by 
communication  links,  especially  when  the  propagation  delay 
r  is  small.  In  other  words,  there  may  not  exist  a  connected 
component  which  has  a  path  with  Ts/t  hops.  Thus,  it  can  be 
expected  that  the  MMD-to-distance  ratio  for  a  network  that  is 
instantaneously  connected  wpp.  is  much  smaller  than  that  for 
one  that  is  not  instantaneously  connected  a.s.  (see  Fig.  5-(c,  d) 
and  Fig.  6-(c,  d)  in  Sec.  V  for  an  illustration). 

V.  Simulation  Results 

We  present  two  sets  of  simulation  results.  One  set  is  to 
show  the  impact  of  connectivity  on  the  scaling  law  of  MMD 
with  respect  to  the  source-destination  distance  (see  Fig.  5  and 
Fig.  6),  and  the  other  set  is  to  show  the  impact  of  the  temporal 
dynamics  of  the  primary  traffic  on  the  MMD  (see  Fig.  7). 
The  density  A 5  of  the  simulated  secondary  network  is  larger 


than  the  critical  density  Ac.  Thus,  the  secondary  network  is 
either  instantaneously  connected  or  only  connected  but  not 
instantaneously  connected,  depending  on  the  density  A pt  of 
the  primary  transmitters.  The  area  of  the  network  is  chosen 
to  be  large  enough  such  that  the  asymptotic  behavior  can  be 
observed.  Without  loss  of  generality,  we  assume  that  the  source 
is  located  at  the  origin.  Each  node  in  the  network  is  a  potential 
destination.  This  allows  us  to  simulate  different  realizations  of 
the  source-destination  pair  using  one  Monte  Carlo  run. 

We  consider  two  mobility  models  for  the  primary  trans¬ 
mitters:  the  random  walk  model  and  the  random  waypoint 
model  [21],  where  the  former  model  has  i.i.d.  increments 
but  the  latter  one  does  not.  For  the  random  walk  model, 
each  primary  transmitter  has  a  probability  p™  of  staying 
at  the  current  position  in  the  next  slot;  otherwise  it  will 
move  to  a  new  position  according  to  a  displacement  vector 
whose  magnitude  is  uniformly  distributed  within  an  interval 
[rm,fif]  and  whose  angle  is  uniformly  distributed  within 
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(c)  instantaneously  connected  (\pt  =  5km  2,  r  =  0.01s)  (d)  not  instantaneously  connected  (Apr  =  30km  2,  r  =  0.01s) 

Fig.  6.  MMD-to-distance  ratio  (in  logarithmic  scale)  vs.  source-destination  distance  for  random  waypoint  model.  Notice  that  the  MMD-to-distance  ratio  is 
obtained  in  one  Monte  Carlo  run.  The  secondary  users  are  distributed  within  a  square  [—5km,  5km]  X  [— 5km,  5km]  with  density  Ay  =  700km-2.  Given  the 
transmission  range  rp  =  50 m  of  the  secondary  users.  Ay  is  larger  than  the  critical  density  Ac(50)  =  570km  2.  Other  simulation  parameters  are:  rj  =  80m, 
Rp  =  50m,  Rj  =  80m,  Tg  =  Is,  pgWp  =  0.05,  vm  =  5m/s,  Vm  =  30m/s. 


[0,  2n).  When  it  reaches  the  simulation  boundary,  it  bounces 
off  the  simulation  border  with  an  angle  determined  by  the 
incoming  direction. 

For  the  random  waypoint  model,  each  primary  transmitter 
chooses  a  random  destination  (not  the  destination  for  its  trans¬ 
mission)  uniformly  distributed  in  the  simulation  area,  which 
determines  its  displacement  direction;  and  then  it  chooses 
a  random  speed  uniformly  distributed  within  an  interval 
to  move  towards  the  destination;  upon  reaching  the 
destination,  it  may  stay  for  a  random  number  of  slots  which  is 
geometrically  distributed  with  parameter  1 — p™p.  The  primary 
receivers  are  uniformly  distributed  within  transmission  range 
Rp  of  their  corresponding  transmitters  in  each  slot. 

Since  it  is  difficult  to  identify  the  path  with  the  MMD  which 
depends  on  the  topology  of  the  secondary  network  and  the 
transmitting  and  receiving  activities  in  the  primary  network 
in  an  intricate  way,  we  obtain  the  MMD  by  considering  a 
flooding  scheme  that  tries  every  possible  path  from  source 


to  destination.  During  flooding,  every  user  that  has  received 
the  message  (including  the  source)  will  transmit  the  message 
to  its  neighbors  (i.e.,  within  its  transmission  range)  when  it 
experiences  a  bidirectional  spectrum  opportunity  with  any  of 
its  neighbors.  The  transmission  attempts  will  not  stop  until  all 
its  neighbors  receive  the  message.  The  time  that  a  user  first 
receives  the  message  during  the  flooding  is  the  MMD  from  the 
source  to  this  user.  To  highlight  the  impact  of  the  waiting  time 
for  spectrum  opportunities  which  is  unique  to  CR  networks,  we 
do  not  consider  the  delay  caused  by  scheduling,  contention,  or 
queuing.  It  can  be  shown  that  this  flooding  scheme  gives  us  the 
MMD.  We  stress  that  flooding  is  used  solely  to  determine  the 
MMD  and  verify  our  scaling  laws;  flooding  is  not  suggested 
as  a  routing  protocol  in  the  secondary  network. 

Fig.  5-(a,  b)  and  Fig.  6-(a,  b)  show  the  MMD-to-distance 
ratio  as  a  function  of  the  source-destination  distance  when  the 
propagation  delay  r  =  0,  where  each  dot  represents  a  realiza¬ 
tion  of  the  destination.  We  see  that  if  the  secondary  network 
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Fig.  7.  The  fractional  difference  (E[MMDi]  —  E[MMD2])/E[MMDi]  expressed  as  a  percentage  vs.  source-destination  distance.  The  average  value  is  taken 
over  1000  Monte  Carlo  runs.  For  E[MMDi],  the  density  A pt(£)  of  the  primary  transmitters  is  fixed  to  be  5km-2.  For  E[MMD2],  A pr(f)  is  uniformly 
distributed  within  [0km-2, 10km-2]  in  each  slot,  and  it  is  i.i.d.  across  slots.  Other  simulation  parameters  are  the  same  as  Fig.  ??. 


is  instantaneously  connected  (Fig.  5-(a)  and  Fig.  6-(a)),  the 
ratio  decreases  rapidly  with  distance  and  can  be  expected  to 
go  to  zero.  On  the  other  hand,  if  the  secondary  network  is  not 
instantaneously  connected  (Fig.  5-(b)  and  Fig.  6-(b)),  the  ratio 
levels  off  as  the  distance  increases  and  will  approach  a  positive 
constant.  Note  that  in  Fig.  5-(a)  and  Fig.  6-(a),  the  MMD-to- 
distance  ratios  of  different  realizations  of  the  destination  are 
grouped  into  several  continuous  curves,  each  associated  with 
a  fixed  MMD.  Since  the  message  is  mainly  delivered  via  the 
ICC  consisting  of  communication  links  when  the  secondary 
network  is  instantaneously  connected,  the  secondary  users  are 
actually  grouped  according  to  the  first  time  that  they  are  in  an 
ICC.  Fig.  5-(a)  and  Fig.  6-(a)  tell  us  that  due  to  the  temporal 
dynamics  of  spectrum  opportunities  caused  by  the  mobility  of 
the  primary  network,  every  node  will  be  part  of  an  ICC  within 
a  few  slots. 

In  Fig.  5-(c,  d)  and  Fig.  6-(c,  d),  we  compare  the  MMD-to- 
distance  ratio  in  a  network  that  is  instantaneously  connected 
and  in  a  network  that  is  not  when  r  is  nonzero  but  small.  The 
four  red  dashed  lines  in  Fig.  5-(c,  d)  and  Fig.  6-(c,  d)  denote 
the  lower  bound  r/rp  imposed  by  the  propagation  delay. 
Although  the  ratio  for  the  network  that  is  instantaneously 
connected  does  not  go  to  zero  due  to  the  nonnegligible 
propagation  delay,  it  is  100  times  smaller  than  the  ratio  for  the 
network  that  is  not.  Notice  that  a  small  group  of  dots  is  located 
at  the  bottom  left  corner  of  Fig.  5-(d).  This  is  because  they  are 
close  to  the  source,  and  their  corresponding  secondary  users 
happen  to  fall  into  the  small  connected  component  formed 
by  communication  links  containing  the  source  in  the  first  few 
slots. 

We  also  compare  the  expected  MMD  (denoted  by 
E[MMDi])  under  a  mobile  primary  network  that  has  fixed 
traffic  load  with  the  one  (denoted  by  E[MMD2])  under  a 
mobile  primary  network  that  has  time-varying  traffic  load8. 
The  fractional  difference  (E[MMDi]  —  E[MMD2])/E[MMDi] 

8For  mobility,  here  we  only  consider  the  random  walk  model. 


is  expressed  as  a  percentage  in  Fig.  7,  where  the  secondary 
network  is  always  instantaneously  connected.  The  fact  that 
the  fractional  difference  is  always  positive  implies  that  the 
introduction  of  another  type  of  temporal  dynamics  reduces 
the  expected  MMD.  Moreover,  when  the  propagation  delay 
t  =  0  (see  Fig.  7-(a)),  the  fractional  difference  is  more  or 
less  constant;  when  r  >  0  (see  Fig.  7-(b)),  it  drops  as  the 
source-destination  distance  increases.  Since  the  percentage  of 
the  secondary  users  in  the  ICC  for  case  2  is  larger  than  the 
one  for  case  1  during  many  slots,  the  waiting  time  of  each 
secondary  user  to  become  part  of  the  ICC,  which  equals  the 
MMD  when  r  =  0,  is  uniformly  decreased  (irrespective  of 
the  distance  from  the  source).  But  when  r  >  0,  the  reduction 
of  the  expected  MMD  in  case  2  is  limited  by  the  positive 
propagation  delay  r. 

VI.  Conclusion  and  Future  Work 

We  say  that  the  secondary  or  cognitive  radio  network  is 
connected  if  the  minimum  multihop  delay  between  an  arbitrary 
source-destination  pair  is  finite  wpp.  We  have  analytically 
characterized  this  connectivity.  The  impact  of  the  primary 
traffic  on  the  connectivity  of  the  secondary  network  has  been 
examined  by  establishing  a  necessary  and  sufficient  condition 
for  connectivity.  Specifically,  depending  on  whether  the  pri¬ 
mary  traffic  has  temporal  dynamics  or  not,  the  connectivity 
of  the  secondary  network  is  equivalent  to  its  topological 
connectivity  which  is  independent  of  the  primary  traffic,  or 
its  instantaneous  connectivity  which  depends  on  the  primary 
traffic.  The  temporal  dynamics  of  the  primary  traffic  can  be 
caused  by  either  mobility  or  changes  in  traffic  load  and  pattern, 
and  it  is  shown  to  significantly  improve  the  connectivity  of 
the  secondary  network  in  the  sense  that  no  matter  how  heavy 
the  primary  traffic  is,  the  secondary  network  is  connected  as 
long  as  its  density  exceeds  the  critical  density  of  homogeneous 
networks. 
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We  have  also  studied  the  impact  of  connectivity  on  the 
multihop  delay.  When  the  propagation  delay  is  negligible,  de¬ 
pending  on  whether  the  secondary  network  is  instantaneously 
connected  with  a  positive  probability  or  not,  the  scaling  of  the 
minimum  multihop  delay  behaves  differently  in  terms  of  the 
scaling  order.  This  scaling  result  is  independent  of  the  random 
positions  of  the  source  and  the  destination,  and  it  only  depends 
on  the  network  parameters  (e.g.,  the  density  of  the  secondary 
users  and  the  traffic  load  of  the  primary  network). 

In  the  above  analysis,  we  have  assumed  a  disk  signal  prop¬ 
agation  model  which  only  incorporates  path-loss.  If  we  take 
into  account  fading  and  shadowing,  then  a  fixed  transmission 
range  does  not  hold,  leading  to  a  random  connection  model 
(RCM)  [3,  Chapter  1],  Since  the  RCM  shares  several  basic 
properties  (e.g.,  the  ergodicity  and  the  existence  of  the  critical 
density)  with  the  Boolean  model  used  in  this  paper,  we  expect 
that  the  results  established  here  can  be  extended  to  the  RCM, 
although  the  derivations  may  become  more  complicated. 

To  highlight  the  impact  of  the  waiting  time  for  spectrum  op¬ 
portunities  on  the  multihop  delay  of  the  secondary  network,  we 
have  not  considered  delay  caused  by  scheduling  or  contention 
among  secondary  users  which  shares  similar  flavor  to  that  in 
conventional  ad  hoc  networks.  The  results  thus  characterize 
the  minimum  delay  and  the  fundamental  performance  limit. 
It  is  our  hope  that  this  paper  will  serve  as  a  starting  point 
to  a  more  complete  characterization  of  multihop  delay  that 
includes  all  these  different  factors. 

Appendix  A:  Basics  of  Ergodic  Theory 

The  study  object  of  ergodic  theory  is  the  so-called  measure¬ 
preserving  (m.p.)  dynamical  system  (d.s.)  (f2,  £F,  n >  T), 
which  consists  of  a  set  f 2,  a  cr-algebra  T  of  measurable  subsets 
of  fl,  a  nonnegative  measure  //  on  (O,  5F),  and  an  invertible 
m.p.  transformation  T  :  fl  — ►  fl  such  that  p(T_1F)  =  n(F) 
V  F  £  9\  A  set  F  £  T  is  said  to  be  T-invariant  if  T~lF  =  F. 

An  m.p. d.s.  (fl,  T,  p,,  T)  is  said  to  be  ergodic  if  for  any 
invariant  set,  either  itself  or  its  complement  has  measure  zero. 
We  use  the  following  fact  frequently  in  the  paper. 

Fact  Al:  [20,  Theorem  2.4.4]  An  m.p. d.s.  (fl,  T,  p,  T) 
is  ergodic,  where  (fl,  5F,  p)  is  a  probability  space,  iff.  V/  £ 
L1(f 1,  T,  p)  (i.e.,  f  is  a  random  variable  with  finite  mean), 
and  u  6  0  we  have 

1  n—  1 

lim  -J2f{Tku)  = 

n — >oo  n  z — ' 
k- 0 

If  T  is  a  shift  transformation  in  the  time  domain,  the  above 
equation  can  be  interpreted  as  the  a.s.  equality  between  the 
time  average  and  the  ensemble  average. 

An  m.p. d.s.  (fl,  5F,  p,  T)  is  said  to  be  mixing  if  V  E,F  £ 
T,  fi(TnE  n  F)  —  h(E)ij.(F)  — ►  0  as  n  — >  oo.  A  mixing 
m.p. d.s  is  ergodic  [20,  Proposition  2.5.1].  Typically  it  is  easier 
to  establish  ergodicity  by  showing  that  the  m.p. d.s.  is  mixing. 

Appendix  B:  Proof  of  Lemma  2 

We  consider  two  cases: 

Case  1:  E[q(f)]  =  E[p(t)]  =  0  but  E[r2(f)]  >  0.  Without 
loss  of  generality,  we  assume  that  E  [m2(i)]  >  0  where 


mx(t )  is  the  x-component  of  the  displacement  m(t)  in  slot 
t.  Consider  two  events  F\  and  F2  which  depend  only  on  the 
points  of  Ilp'r(O)  within  a  box  Bm  centered  at  the  origin 
with  side  length  m.  Let  Gn  denote  the  event  that  all  the 
points  of  IIp7’(0)  within  Bm  are  outside  Bm  in  slot  n,  and 
Hk  the  event  that  there  are  at  most  K  points  of  IIpp(O) 
within  Bm.  Let  X(n)  be  the  .'/.'-component  of  the  cumulative 
displacement  associated  with  a  point  from  t  =  0  to  t  =  n,  i.e., 
X(n )  =  J2t= omx{t)-  Then 

Pr{Gn\HK} 

K 

>  (Pr{ | JC (?x  —  1)|  >  m})k  Pr{#points  £  Bm  =  k\HK} 

k—0 

>  (Pr{|X(n  —  1)|  >  m})K  . 

Here  we  have  used  the  i.i.d.  property  of  X(t)  across  points. 

It  follows  from  Lemma  B1  in  [22,  Appendix  B]  that 

lim  Pr{|X(n  —  1)|  >  m}  =  1. 

n — >oo 

Thus, 

lim  Pv{Gn\HK}  =  1.  (Bl) 

n — xx) 

If  Gn  occurs  for  some  n,  then  obviously  F\  is  independent 
of  TnF2,  i.e.,  Pr{F!  n  TnF2\Gn  n  HK}  =  Pr{F1|Gn  n 
HK}Pr{TnF2\Gn  n  Hk}.  Now  since 

Pr { F\  n  Tn F2\Hk} 

=  Pr{Fi  n  TnF2  | Gn  n  HK}Pr{Gn\HK) 

+Pr{Fi  n  TnF2\G ~  n  HK}Pr(G^\HK}, 

we  have  that 

P*{Fi\Gn  n  HK}Px{TnF2\Gn  n  HK}Pr{Gn\HK} 

<  Pr{Fi  fl  Tn F2\H k} 

<  PxiF^Gn  n  HK}Px{TnF2\Gn  n  HK}Px{Gn\HK} 
+Px(G^\Hk}. 

Thus, 

lim  Pr{Fi  n  TnF2\HK] 

n — >oo 

=  lim  Pr{Fi|Gn  n  HK}Px{TnF2\Gn  n  HK}Px{Gn\HK} 

n — >00 

=  lim  Pr{Fi |G„  n  HK}  lim  Px{TnF2\Gn  n  HK}. 

n — >oo  n — >oo 

Eqn.  (Bl)  and  the  temporal  stationarity  of  { 1 1 i>t(F)}  yield 
lim  Px{Fi\Gn  PI  Hk}  =  Pr{Fi|ff*}, 

n — >oo 

lim  Px{Tn F2\Gn  fl  Hk}  =  lim  Px{TnF2\HK} 

n — >oo  n — >oo 

=  Px{F2\Hk}. 

We  thus  have  that 

lim  Pr{Fi  fl  TuF2\Hk}  =  Pr{Fi \HK}Px{F2\HK}- 

n — »-oo 

Since  lim  Px{Hk}  =  1,  as  K  — >  oo, 

K  — >oo 

lim  Pr{Fi  n  TnF2}  =  Pr{Pi}Pr{F2}. 

n — >oo 

Since  any  two  arbitrary  events  F\  and  F2  can  be  approxi¬ 
mated  by  two  sequences  of  events  {F™}  and  {F™}  which 


/ d/i  a.s. 
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depend  only  on  the  realization  of  IIpt(O)  inside  Bm,  the 
conclusion  follows  from  Lemma  Cl  in  [22,  Appendix  C]. 

Case  2:  E[g(£)]  >  0.  Consider  two  events  F\  and  F2  which 
depend  only  on  the  points9  of  {IIpt(£)  :  —  T  <  t  <  T} 

within  a  box  Bm  centered  at  the  origin  with  side  length  rri. 
Let  Gn  denote  the  event  that  all  the  points  which  have  visited 
Bm  during  —T  <  t  <  T  do  not  transmit  in  slot  n  —  T,  and 
Hk  the  event  that  there  are  at  most  K  such  points.  Fixing  a 
realization  of  {g(f}},  we  have  that 


> 


> 


Pr{Gn\HK} 

K 

E 

k= 0 


i-T- 1 


1  -  n  (1  -?00) 


i=T 


Pr{#points  E  Bn 


i-T-l 


l  n  (l-g(*)) 


i=T 


K 


k\HK} 


Since  {g(f)}  is  ergodic,  based  on  Fact  Al,  we  have  that 


n—  1 


iim  -  V  q(t)  =  E[g(f)]  a.s., 

n— >  00  77,  z ' 
t= 0 

implying  that  lim  q{t)  =  00  a.s.  It  follows  that 


n-T-l 


(3).  It  follows  from  the  arguments  similar  to  those  in  showing 
that  tsw  <  00  a.s.  in  the  proof  of  Lemma  1,  that  to  <  00  a.s. 

Given  C(to),  we  define  user  wv  as  the  user  in  C'(fo)  that 
is  closest  to  v ,  i.e., 

wv  =  argmin  d(wi,  v). 
wieC(t0) 

Notice  that  if  v  G  C(to),  then  wv  =  v. 

As  illustrated  in  Fig.  4,  the  constructed  path  Lc  passes 
through  wv,  and  the  multihop  delay  tc(fi,v)  along  the  path 
Lc  can  be  expressed  as 

tC(n,  v)  =  t0+  t(fi,  wv)  +  t(wv,u)  =t0  +  t(wv ,  v), 

where  t(wv,v)  is  the  MMD  from  wu  to  v.  In  the  last  step, 
we  have  used  t(ft,wv)  =  0,  which  is  due  to  the  fact  that  /i, 
wv  6  C'(fo)  and  t  =  0.  Next  we  prove  the  following  lemma. 

Lemma  Cl:  t{wv,v)  is  finite  a.s. 

Proof  Sketch:  We  first  show  that  d{wv,v)  <  00  a.s.  by 
using  the  ergodicity  of  the  network  model,  and  then  obtain 
an  upper  bound  on  the  multihop  delay  tL  (wv ,  v)  along  the 
shortest  path11  L(wv,v)  from  wv  to  v.  Since  t{wv,v)  < 
tL(wv,  v),  the  a.s.  finiteness  of  t(wv,  v)  follows  from  that  of 
the  upper  bound  on  tL(wv,  v).  The  proof  here  is  inspired  by 
the  proof  of  Lemma  9  in  [13],  but  with  a  much  simpler  proof 
of  d(wv,  v)  <  00.  For  details,  see  [15,  Appendix  B].  ■ 


lim  TT  (1  —  q{i))  =  0  a.s. 

n — >00 

i=T 


Appendix  D:  Proof  of  T2.2 


Thus, 

lim  Pr {Gn\HK}  =  1, 

n — ►oo 

and  the  rest  of  the  proof  follows  along  the  same  line  of  the 
one  of  Case  1. 


Let  t'(n,  v)  be  the  MMD  from  /1  to  v  when  Xprit)  =  X'PT 
V  t.  Then  based  on  a  coupling  argument  [3,  Chapter  2],  we 
have  E[f(/x,  v)\  >  F.[t'(n,v)\.  It  suffices  to  show  that 


lim  inf 

— ►00 


E[£'(/t,^)] 

d(p,  v) 


Appendix  C:  Proof  of  T2. 1 


which  is  shown  in  [15,  Lemma  4], 


We  use  the  ICC  consisting  of  communication  links  in 
Qh (As,  Apt(£o)>  fo)  during  some  slot  to  to  construct  a  path 
Lc  from  //  to  v  such  that  the  multihop  delay  along  this  path 
is  independent  of  the  distance  d(/r,  v)  (see  Fig.  4).  Then  we 
analyze  the  multihop  delay  tc(n,  v)  along  Lc- 

Assume  that  ft  starts  trying  to  send  the  message  at  time  t  = 
0.  Let  C(t)  be  the  ICC  in  Qh ( As ,  Xpr{t),t)  if  it  exists 10,  and 
to  the  first  slot  such  that  ft  G  C(to).  Let  (Dp,  5Fp,  Pp)  be  the 
probability  space  of  {Ilpp(f),  IIpp(£)  :  t  G  Z},  and  define  a 
shift  transformation  T  similarly  to  (4).  Given  that  E[q(f)]  >  0 
or  E[r2(f)]  >  0,  the  ergodicity  of  the  m.p.  dynamical  system 
{DpjfFp,  Pp,T}  follows  along  the  same  line  of  the  proof 
of  Lemma  2.  Let  Ft  denote  the  event  that  ft  G  C(t).  Since 
Pr{App(f)  <  ApT(As)}  >  0  implies  that  3  X'PT  <  XPT(X s) 
such  that  Pr{App(£)  <  X'PT}  >  0,  we  have  that 

Pr{F0}  >  Pr{App(£)  <  X'PT}0(Xs ,  X'PT)  >  0, 

where  6(Xs,X'PT)  is  the  probability  that  an  arbitrary  sec¬ 
ondary  users  belongs  to  an  ICC  in  Qh(Xs,  X'PT,t)  given  by 

^  Si  nee  the  set  of  the  primary  transmitters  may  change  in  every  slot,  it  is 
not  enough  to  only  consider  the  points  of  npp(0). 

nSincc  App(t)  is  time-varying,  it  is  possible  that  C(t)  does  not  exist  in 
some  slots. 
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